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Abstract 

We examine relations among the parameters characterizing the phenomenologi- 
cal equation of state (EOS) of nearly symmetric, uniform nuclear matter near the 
saturation density by comparing macroscopic calculations of radii and masses of 
stable nuclei with the experimental data. The EOS parameters of interest here are 
the symmetry energy So, the symmetry energy density-derivative coefficient L and 
the incompressibility Kq at the normal nuclear density. We find a constraint on the 
relation between K$ and L from the empirically allowed values of the slope of the 
saturation line (the line joining the saturation points of nuclear matter at finite 
neutron excess), together with a strong correlation between So and L. In the light 
of the uncertainties in the values of Kq and L, we macroscopically calculate radii 
of unstable nuclei as expected to be produced in future facilities. We find that the 
matter radii depend strongly on L while being almost independent of K$, a feature 
that will help to determine the L value via systematic measurements of nuclear size. 
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1 Introduction 



Saturation of density and binding energy, one of the fundamental properties 
of atomic nuclei, underlies a liquid-drop approach described by a Weizsacker- 
Bethe mass formula [1] 

— Eb = E YO \ + E sym + E SUT { + Ecou\, (1) 
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where Eb is the nuclear binding energy, E vo \ is the volume energy, E sym is the 
symmetry energy, -E sur f is the surface energy, and -Ecoui is the Coulomb energy. 
The sum E vo \ + E syia corresponds to the saturation energy of uniform nuclear 
matter. Since matter in nuclei is a strongly interacting system, it remains a 
challenging theoretical problem to understand the nuclear matter equation 
of state (EOS) through microscopic calculations that utilize a model of the 
nuclear force duly incorporating low-energy two-nucleon scattering data and 
properties of light nuclei [2]. Furthermore, it is not straightforward to empir- 
ically clarify the EOS, although constraints on the EOS from nuclear masses 
and radii (e.g., Refs. [3-5]), observables in heavy-ion collision experiments per- 
formed at intermediate and relativistic energies (e.g., Refs. [6-9]), the isoscalar 
giant monopole resonance in nuclei (e.g., Ref. [10]) and even X-ray observa- 
tions of isolated neutron stars [11] are available. In this work we will consider 
such constraints from nuclear masses and radii. 



The EOS of bulk nuclear matter is a function of nucleon density n and pro- 
ton fraction x, which are related to the neutron and proton number densities 
n n and n p as n n = n(l — x) and n p = nx. We may generally express the energy 
per nucleon near the saturation point of symmetric nuclear matter as [12] 



w = w + 



K 
18n 2 n 



(n - n ) 2 + 
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3n 



a 



(2) 



Here w , n and K are the saturation energy, the saturation density and 
the incompressibility of symmetric nuclear matter, So is the symmetry energy 
S(n) at n = no, L = 3no(dS / dn) n=no is the density symmetry coefficient, and 
a = 1 — 2x is the neutron excess. As the neutron excess increases from zero, 
the saturation point moves in the density versus energy plane. This movement 
is determined mainly by the parameters L and So associated with the density- 
dependent symmetry energy S(n) [3]. Up to second order in a, the saturation 
energy w s and density n s are given by 

w s = w + S a 2 (3) 



and 
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3n L 
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The slope, y, of the saturation line near a — (x — 1/2) is thus expressed as 

KqSq 

y = -3n^L- (5) 
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Derivation of L and So from nuclear observables is generally obscured by 
the interfacial and electrostatic properties. Among the observables, the masses 
and root-mean-square radii of nuclei, which are controlled mainly by the bulk 
properties, are expected to be good tracers of L and S - This expectation 
was stressed by an earlier investigation [3] based on two macroscopic nuclear 
models. Such macroscopic models are reliable in a range of neutron excess, 
a < 0.3, and mass number, A > 50. In this range the neutron separation 
energy that can be evaluated from a Weizsacker-Bethe mass formula is greater 
than 2 MeV, allowing us to preclude the possibility of neutron halo formation 
expected at large neutron excess (or small separation energy). In constraining 
the EOS from masses and radii of nuclei of neutron excess a < 0.3 and mass 
number A > 50 within the framework of the macroscopic models, systematic 
study allowing for uncertainties in the EOS parameters is indispensable. 

In this paper we thus explore a systematic way of extracting L and So from 
empirical masses and radii of nuclei, together with the parameters, n , w and 
K , characterizing the saturation of symmetric nuclear matter. We first set 
an expression for the energy of uniform nuclear matter, which reduces to the 
phenomenological form (2) in the limits of n — > no and a — > {x — > 1/2). Using 
this energy expression within a simplified version of the extended Thomas- 
Fermi approximation, which permits us to determine the macroscopic features 
of the nuclear ground state, we calculate charges, charge radii and masses of 
/5-stable nuclei. Comparing these calculations with empirical values allows us 
to derive the optimal parameter set for various values of the slope y and the 
incompressibility Kq. We thus find a strong correlation between L and So- The 
next step is to calculate root-mean-square charge and matter radii of more 
neutron-rich nuclei that are expected to be produced in future radioactive ion 
beam facilities. The results suggest that the density symmetry coefficient L 
may be constrained from possible systematic data on the matter radii in a way 
nearly independent of the poorly known Kq, while the slope y being deducible 
as a function of K . We finally discuss an empirically allowed region in the 
space of the parameters characterizing the EOS (2). 

In Section 2 we construct a macroscopic model of nuclei. Optimization as- 
sociated with fitting to empirical data for nuclei on the smoothed /5-stability 
line is illustrated in Section 3. In Section 4 we calculate matter and charge 
radii of unstable nuclei. Our conclusions are presented in Section 5. Numerical 
tables are given in the Appendix. 
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2 Macroscopic description of nuclei 



In constructing a macroscopic nuclear model, we begin with the expression 
for the bulk energy per nucleon [13], 

Ofc2/o 2N2/3 

w = \ 1 (nf 3 + nf) + (1 - a 2 )v s (n)/n + a 2 v n (n)/n, (6) 



where 

,3 



v s = a x n + — (7) 

1 + a 3 n 



and 
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are the potential energy densities for symmetric nuclear matter and pure neu- 
tron matter, and m n is the neutron mass. (Replacement of the proton mass 
m p by m n in the proton kinetic energy makes only a negligible difference.) For 
the later purpose of roughly describing the nucleon distribution in a nucleus, 
we incorporate into the potential energy densities (7) and (8) a low density 
behaviour oc n 2 as expected from a contact two-nucleon interaction. Note, 
however, that we will focus on the EOS of nearly symmetric nuclear matter 
near the saturation density. We will thus determine the parameters included 
in Eqs. (7) and (8) in such a way that they reproduce data on radii and masses 
of stable nuclei. 

In the limit of n — > no and a — > (x — > 1/2), expression (6) reduces to the 
usual form (2) according to 

1 /3vr 2 \ 2/3 £ 2 2/3 ( b 2 a 2 \ 2 

Q \ 2 J m n \1 + b 3 n 1 + a 3 n J 
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Equation (13) comes from the fact that the density derivative of the energy 
per nucleon at n = n and a = (x — 1/2) vanishes due to the saturation. 
Note that the five relations (9)— (13) are not sufficient to determine the six 
parameters ai,...,6 3 , and that the parameter b 3 , which controls the EOS 
of matter at large neutron excess and high density, has little effect on the 
saturation properties of nearly symmetric nuclear matter. We will thus set 63 
as a typical value 1.58632 fm 3 , which was obtained by one of the authors [13], 
and determine the other parameters from the empirical radii and masses of 
stable nuclei. 

We now proceed to describe a spherical nucleus of proton number Z and 
mass number A within the framework of a simplified version of the extended 
Thomas- Fermi theory [13]. We first write the total energy of a nucleus as a 
function of the density distributions n n (r) and n p (r) according to 

E = E b + E g + E c + Nm n + Zm p , (14) 

where 

E b = J d 3 rn(r)w(n n (r),n p (r)) (15) 

is the bulk energy, 

E g = F j d 3 r\Vn(r)\ 2 (16) 

is the gradient energy with an adjustable constant F , 

Ec= e lj^J^-f^l (17) 



is the Coulomb energy, and N = A — Z is the neutron number. Here we 
ignore shell and pairing effects. We also neglect the contribution to E g from 
the gradient of the proton fraction x [see Eq. (16)]; this contribution makes 
only a little difference even in the description of extremely neutron-rich nuclei, 
as clarified in the context of neutron star matter [13]. 
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The energy (14), once optimized, can be mapped onto a Weizsacker-Bethe 
mass formula of the form (1) via 



E SU]: { = E g + (E b — E vo \) + (E c — E C oul), 



(18) 



where E b — E vo i denotes the inhomogeneity contribution to the bulk energy E b , 
and E c — E Co u\ (E Co u\ = SZ 2 e 2 /5R with the liquid-drop radius R) denotes that 
to the Coulomb energy Ec- These inhomogeneity contributions arise from the 
fact that matter in a nucleus is compressible. We remark that in equilibrium 
with respect to nuclear size, E g = Ec holds [13]. Combining this relation 
with the well-known equilibrium condition for the liquid-drop size, E sur f = 
2£'coui, an d with -Ecoui — Ec, we find a simple relation, E g ~ E SUT f/2, for an 
equilibrium nuclide. 

For the present purpose of examining the macroscopic properties of nuclei 
such as masses and radii, it is sufficient to characterize the neutron and pro- 
ton distributions for each nucleus by the central densities, radii and surface 
diffuseness different between neutrons and protons, as in Ref. [13]. We thus 
assume the nucleon distributions n«(r) (i = n,p), where r is the distance from 
the center of the nucleus, as 



Here Ri roughly represents the nucleon radius, the relative surface diffuse- 
ness, and nf the central number density. The proton distribution of the form 
(19) can fairly well reproduce the experimental data for stable nuclei such as 
90 Zr and 208 Pb, as we shall see in the next section. 



3 Optimization 

For fixed mass number A, we then optimize the energy (14) with respect 
to the parameters Ri, t« and n- n . (In a general Thomas- Fermi approach, such 
optimization is carried out without assuming a particular form of the distri- 
butions.) The resultant optimal values of charge number, nuclear mass and 
root-mean-square charge radius 



rii(r) = < 




(19) 
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Fig. 1. Various optimal relations among the parameters Sq, uq, wo, L and Kq 
characterizing the EOS of nearly symmetric nuclear matter. In addition to the 
present results (crosses), the Skyrme-Hartree-Fock predictions [dots except for SIII 
(square)] and the TM1 prediction (triangle) are plotted. In (c), the thin lines are 
lines of constant y. 



where 



p c (r) = (tt 1 / 2 ^)- 3 / rfVexp (-|r - rf/a 2 p ) n p (r>) 



(21) 



with a p = 0.65 fm is the charge distribution folded with the proton form factor 
[14], are functions of ai-63 and Fq. These optimal values are in turn compared 
with the empirical values for nuclei on the smoothed /^-stability line ranging 
25 < A < 245 (see Table A.l in Ref. [13], which is based on Refs. [15,16]). For 
fixed slope y and incompressibility K , such a comparison can be made by a 
usual least squares fitting, which gives rise to an optimal set of the parameters 
ai-63 and F . Here, we set y and K as -1800 MeV fm 3 < y < -200 MeV fm 3 
and 180 MeV < K < 360 MeV; the numerical results for n , u>o, Sq, L and F 
are tabulated in the Appendix. When the slope is very gentle (0 > y > —200 
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Fig. 2. The energy per nucleon as a function of nucleon density. The dotted lines 
denote the saturation lines. 

MeV fm 3 ), the optimal set is unavailable for large K ; when the slope is very 
steep, the optimal parameters converge on the values close to those obtained 
for y = —1800 MeV fm 3 . Nuclear masses that can be calculated from the 
optimal parameter sets agree well with the experimental data of 1962 nuclides 
(A > 2) [17]; for all the combinations of y and Kq, the root-mean-square 
deviations of the masses are ~ 3-5 MeV, which are roughly as large as those 
obtained from a Weizsacker-Bethe type formula. We likewise evaluated the 
root-mean-square charge radii of various stable nuclei; the root-mean-square 
deviations from the experimental data of 92 nuclides (A > 50) [16] are about 
0.06 fm. 

The optimal results for So, L, n and wo are plotted in Fig. 1. From Fig. 
la, we obtain a relation nearly independent of K : 

S « B + CL, (22) 
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where B m 28 MeV and C m 0.075. We find from Figs, lb and Id that the 
saturation energy and density of symmetric nuclear matter always take on a 
value of — 16.0±0.5 MeV and 0.155±0.015 fm~ 3 . Several data in Fig. Id having 
n smaller than the standard values correspond to the case of y — —200 MeV 
fm 3 , where the fitting is no longer effective. In Fig. lc, we see a band on which 
the optimal values of L and Kq are scattered; in this band L increases with 
increasing y for fixed Kq. For comparison we also plot the predictions from 
various Skyrme-Hartree-Fock schemes (references in Ref. [4]; Refs. [18-20]) 
and a relativistic mean field model (TM1 in Ref. [21]); the values for n , w , 
K , S , L and y are tabulated in Table I. These predictions are distributed 
over the band, among which only two (TM1 and SIII) were considered in the 
previous analysis [3]. We remark in passing that the optimal values of Fq are 
confined to 66 ± 6 MeV fm 5 , consistent with the result of Ref. [13]. 

In Fig. 2 we display the EOS (6) for various sets of y and K . Whereas y 
affects the slope of the saturation line, Kq controls tiq as well as the curvature 
of the line of constant a. The neutron and proton distributions in 208 Pb and 
90 Zr modelled via Eq. (19) are plotted in Fig. 3. The question we consider in 
the next section is how such differences in the saturation properties as shown 
in Fig. 2 affect matter and charge radii of unstable nuclei that can be evaluated 
from the distributions of the form illustrated in Fig. 3. 



4 Matter and charge radii 

For neutron-rich nuclides we now obtain the root-mean-square charge radii 
R c and matter radii R m , defined as 

1/2 

(23) 
where 

Pm (r) = (ir^X)- 3 ! dVexp (-|r - rf/aj) n(r') (24) 

is the matter distribution folded with the proton charge form factor equally for 
neutrons and protons. We evaluated the radii R c and R m of Ni and Sn isotopes 
for combinations of y = -220, -350, -1800 MeV fm 3 and K = 180, 230, 360 
MeV. The results are shown at neutron excess of up to a ~ 0.4 in Fig. 4. We 
find from the upper panels of Fig. 4 (see also the upper panels of Fig. 5 for 
further clarity) that at a ~ 0.3 (x ~ 0.35) a difference of order 0.05-0.1 fm 
occurs in the matter radii due to variation in the slope y. This is because as 
the slope becomes gentler, the saturation density difference n — n s becomes 
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Fig. 3. The neutron and proton density distributions in 208 Pb and 90 Zr. 
Experimental data on the charge distributions (crosses) are taken from Ref. [16]. 



larger [see Eq. (4)]. We also find that the charge radii depend only weakly 
on y, a feature consistent with the prediction made in Ref. [3]. We can thus 
expect that forthcoming empirical data on matter radii of unstable neutron- 
rich nuclei with accuracy down to order ±0.01 fm will at least answer the 
question of whether the slope y is steep or gentle. 

Several remarks, however, are needed here. First, it is to be noted that the 
experimental matter radii that can be derived from measurements of interac- 
tion cross sections [22] and elastic scattering of protons and alpha particles 
[23] depend strongly on treatment of the optical potential, in contrast to the 
charge radii that can be determined from elastic electron scattering [16,24], 
muonic X-ray experiments [16,24] and isotope-shift measurements [25,24]. This 
dependence contributes to intrinsic uncertainties in the derived matter radii, 
which are typically of order or greater than ±0.05 fm [23,26]. Second, we recall 
that the present calculations of the radii R m and R c ignore the tails of the 
nucleon distributions arising from quantum-mechanical effects; the absence 
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Fig. 4. The root-mean-square charge and matter radii of Ni and Sn isotopes for 
combinations of y = -220, -350, -1800 MeV fm 3 and K = 180, 230, 360 MeV. 
Experimental data on the root-mean-square charge radii (dots) and matter radii 
(crosses) are taken from Refs. [16] and [23], respectively. 

of such tails tends to reduce the radii. Third, there are uncertainties in the 
calculated radii due to the absence of shell and pairing effects in the present 
macroscopic models. These models, which are fitted to the charge radii of nu- 
clei on the smoothed f3 stability line, provide the proton-closed-shell nuclei, Ni 
and Sn, with larger charge radii than the empirical values, as shown in Fig. 4. 
Fourth, we can see from the lower panels of Fig. 4 that with Kq increased and 
y fixed, the matter radii R m increase, a feature that prevents a clear derivation 
of y from R m . 

Such Kq dependence of the matter radii R m is due mainly to the K de- 
pendence of the saturation density, n s , given by Eq. (4). For fixed y, the K 
dependence of n s is dominated by n , which increases with decreasing K as 
in Fig. Id. Generally, this increase in hq tends to increase the inner nucleon 
densities (see the lower panels of Fig. 3) and hence reduce the radii R rn . Note, 
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Fig. 5. Differences of the root-mean-square matter radii of Ni and Sn isotopes, 
AR m , from those calculated for y = —350 MeV fm 3 and Kq = 230 MeV (upper 
panels) and for L = 50 MeV and Kq = 230 MeV (lower panels). The thin curves in 
the lower panels are from formulas (26) and (27). 

however, that there is an opposite effect of Kq on R m . This effect comes from 
the fact that a significant part of the nucleons are present in the deepest region 
of the nuclear surface where r 2 n n (r) and r 2 n p (r) are peaked. In this region, 
the nucleon densities begin to drop in such a way that with decreasing K , the 
surface diffuses further away. This diffuseness, which can also be seen from the 
lower panels of Fig. 3, is consistent with the results from microscopic nuclear 
models [27]. 

A cancellation between those counteracting effects is favoured for the pur- 
pose of deriving information about the saturation properties in a way inde- 
pendent of Kq. It turns out that better cancellation can be achieved if we 
calculate, as in the lower panels of Fig. 6, the matter radii of Ni and Sn iso- 
topes for the EOS parameters optimized under fixed L rather than y. The 
charge radii likewise calculated are also nearly independent of K . On the 
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Fig. 6. Same as Fig. 4 for combinations of L = 5, 50, 80 MeV and Kq 
360 MeV. 



180, 230, 



other hand, under fixed Kq, the dependence of the matter and charge radii 
on L, as depicted in the upper panels of Fig. 6, is the same as that on y. 
These results from Fig. 6 may open an opportunity to empirically determine 
the density symmetry coefficient L from the isotopic dependence of matter 
radii. 



For the purpose of parametrizing the matter and charge radii as functions of 
A, a and L, it is useful to first obtain fitting formulas for the root-mean-square 
neutron and proton radii, defined as 



Rn 



N' 1 j d 3 rr 2 n n (r) 



1/2 



Rr. 



Z- 1 J d 3 rr 2 n p {r) 



1/2 



(25) 



To be fitted to are the radii R n and R p of Ni, Sn and Pb isotopes ranging 
< a < 0.3 (0.35 < x < 0.5) as calculated from the present macroscopic 
model for the EOS parameter sets tabulated in Table A.l. In Fig. 7 the results 
for 56 Ni, 80 Ni, 116 Sn and 142 Ni are plotted as functions of L. We observe from 
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Fig. 7. The root-mean-square neutron (dots) and proton (crosses) radii of 58 Ni, 
80 Ni, 116 Sn and 142 Sn for the EOS parameter sets tabulated in Table A.l. The 
solid and dotted lines are from formulas (27) and (26), respectively. 

this figure that no uncertainties in R, p more than ±0.03 fm arise from the 
various values of K and L, while the K dependence of R n can be ignored as 
compared with the L dependence, which is stronger at larger neutron excess. 
These features, consistent with the dependence of R m and R c on L and K$ as 
mentioned above, lead us to the formulas, 



R p = CxA 1 ^ + c 2 + c 3 (a - a ) 2 , 



(26) 



with ci = 0.914961 fm, c 2 = -0.102372 fm, c 3 = 0.388905 fm and a = 
0.879704, and 



R n = c 4 A 1/3 (l + c 5 La 2 + c 6 L 2 a 4 ) + c 7 + c s a, 



(27) 



with c 4 = 0.880489 fm, c 5 = 0.00635080 MeV" 1 , c 6 = -0.000172275 MeV" 2 , 
c 7 = 0.301616 fm and c 8 = 0.193326 fm. Expressions (26) and (27) reproduce 
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the original values with the root-mean-square deviations of 0.010 fm and 0.013 
fm, respectively, and roughly obey a usual A 1 / 3 law for fixed L and a. Finally, 
the parametrization of the matter and charge radii can be constructed from 
Eqs. (26) and (27) via the relations R 2 m = (Z/A)R 2 p + (N/A)R 2 n + 3a 2 p /2 and 
R\ = R 2 p + 3a 2 p /2 that can be derived from Eqs. (20), (23) and (25). 

The lower panels of Fig. 5 depict the matter radii of Ni and Sn isotopes 
calculated for various combinations of L and K relative to those calculated 
for L = 50 MeV and K = 230 MeV. As can be seen from these panels, K - 
induced uncertainties in the matter radii at constant L are limited to the order 
of ±0.015 fm in the neutron-rich side, small enough for the L dependence of 
the matter radii to be seen clearly. This may play a role in determining L from 
measurements of the matter radii. Note, however, that possible data on matter 
radii for isotopes of a specific element would be insufficient to determine the 
density symmetry coefficient L, as inferred from the differences between the 
empirical and calculated charge radii shown for Ni and Sn in Fig. 6. In view of 
such differences as well as intrinsic uncertainties in the matter radii that can 
be deduced from proton and alpha-particle elastic scattering data, systematics 
of the isotopic dependence of matter radii would be required. 

It is interesting to note that the difference between the neutron and proton 
radii increases with increasing L in a way dependent on neutron excess but 
almost independent of K (see Fig. 7). This leads us to study the isotopic 
dependence of neutron skin thickness in order to derive the density symme- 
try coefficient L. The previous investigation [3] has already suggested the 
possibility that the EOS of asymmetric nuclear matter will be probed by fu- 
ture detection of neutron skin thickness in unstable neutron-enriched nuclei. 
However, the neutron skin structure depends strongly on the EOS of nuclear 
matter at large neutron excess and low density [4], which is empirically hard 
to determine. Moreover, the fact [28] that the neutron skin thickness is ther- 
modynamically relevant to the dependence of the surface tension on neutron 
excess complicates an extraction of the bulk properties of asymmetric nuclear 
matter; the connection between the surface and bulk properties remains to 
be clarified. We will elsewhere reinvestigate how the neutron skin thickness is 
phenomenologically related to the EOS of nearly symmetric nuclear matter 
[29]. 



5 Conclusion 

In this paper we have derived the relations between the parameters charac- 
terizing the EOS of nearly symmetric nuclear matter from experimental data 
on radii and masses of stable nuclei. We have found the linear relation (22) 
between the parameters L and S associated with the symmetry energy. It is 
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interesting to regard the band structure on the K versus L plane (see Fig. lc) 
as an empirically allowed region of K and L. Future systematic measurements 
of the root-mean-square matter radii of unstable neutron-rich nuclei may nar- 
row the allowed region in such a way as to fix L almost independently of K . 
Once L is determined, the slope y would be given as a function of Kq. It is in- 
structive to note that the incompressibility was estimated to be K = 231 ± 5 
MeV from the observations of the giant monopole resonances [10]. This esti- 
mate can stringently limit the allowed region, although it is model-dependent 
in the sense that it involves microscopic calculations of the resonance energies 
from a specific model of the effective nucleon-nucleon interactions [30]. 
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Appendix: Evaluation of the parameters in Eq. (2) 



In this appendix we tabulate (Table A.l) the optimal values of various EOS 
parameters in Eq. (2). 
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Table A.l 

Optimal values of the parameters no, wq, So, L and Fo for 
fixed y and Kq 
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Table A.l — continued 
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Table A.l — continued 
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Table A.l — continued 
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31, 


.627 


58, 


.000 


64. 


364 


400 


340 


0. 


14938 


16. 


.164 


31. 


.976 


60, 


.648 


67. 


.165 


400 


350 


0. 


14854 


16 


.071 


31, 


.895 


62. 


.628 


63. 


.808 


400 


360 


0. 


14701 


16. 


.095 


32 


.082 


65. 


.471 


65. 


.588 



22 



Table A.l — continued 



-y 


K 




n 


- 


w 


So 


L 


^0 


500 


180 


0. 


16895 


16 


.213 


29. 


.750 


21 


,131 


70. 


.289 


500 


190 


0. 


16674 


16. 


.194 


29 


,786 


22, 


.628 


69, 


.745 


500 


200 


0. 


16506 


16. 


.201 


29 


,868 


24, 


,127 


69, 


.843 


500 


210 


0. 


16324 


16 


.179 


29. 


.901 


25 


.644 


69. 


044 


500 


220 


0. 


16136 


16 


.149 


29 


,914 


27, 


,191 


68. 


.228 


500 


230 


0. 


15991 


16. 


,141 


29 


,972 


28, 


,739 


67. 


.851 


500 


240 


0. 


15856 


16 


.130 


30 


.028 


30 


.302 


67. 


.377 


500 


250 


0. 


15731 


16. 


.120 


30 


,085 


31 


,875 


66. 


940 


500 


260 


0. 


15620 


16 


.112 


30 


.152 


33 


.459 


66. 


.504 


500 


270 


0. 


15513 


16 


.104 


30 


,216 


35 


.060 


66. 


079 


500 


280 


0. 


15421 


16. 


.098 


30 


,288 


36 


,663 


65. 


,685 


500 


290 


0. 


15330 


16 


.092 


30 


.359 


38. 


.286 


65. 


293 


500 


300 


0. 


15243 


16 


.086 


30 


.437 


39 


.934 


64. 


956 


500 


310 


0. 


15165 


16 


.080 


30 


.511 


41. 


,581 


64. 


,550 


500 


320 


0. 


15091 


16 


.076 


30 


,594 


43. 


,249 


64. 


,240 


500 


330 


0. 


15022 


16 


.073 


30 


.678 


44. 


.930 


63. 


950 


500 


340 


0. 


14957 


16. 


.069 


30 


,765 


46. 


,624 


63 


,597 


500 


350 


0. 


14892 


16 


.063 


30. 


.852 


48. 


,341 


63 


250 


500 


360 


0. 


14912 


16 


.222 


31, 


,314 


50. 


.399 


68. 


034 


600 


180 


0. 


16890 


16 


.208 


29 


,461 


17. 


,443 


70. 


,159 


600 


190 


0. 


16669 


16 


.188 


29 


,475 


18. 


.664 


69. 


602 


600 


200 


0. 


16470 


16 


.172 


29 


.496 


19 


,899 


69. 


132 


600 


210 


0. 


16393 


16 


.245 


29 


,692 


21 


.132 


70. 


991 


600 


220 


0. 


16137 


16 


.148 


29. 


.560 


22. 


.389 


68. 


204 


600 


230 


0. 


15990 


16 


.135 


29. 


.596 


23 


.651 


67. 


,684 


600 


240 


0. 


15856 


16. 


,123 


29 


,634 


24. 


.920 


67. 


,189 


600 


250 


0. 


15735 


16 


.114 


29 


.675 


26 


,194 


66. 


,714 


600 


260 


0. 


15626 


16 


,108 


29 


,724 


27. 


.476 


66. 


,344 


600 


270 


0. 


15519 


16 


,097 


29 


,766 


28. 


.770 


65. 


,856 
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Table A.l — continued 



-y 


K 




no 


- 


Wo 


S 




L 


Fo 


600 


280 


0. 


15422 


16 


.091 


29 


.818 


30 


.077 


65 


.479 


600 


290 


0. 


15336 


16. 


.085 


29 


.871 


31 


.381 


65. 


,049 


600 


300 


0. 


15254 


16 


.081 


29. 


.930 


32 


.701 


64 


.714 


600 


310 


0. 


15176 


16. 


.075 


29 


.988 


34. 


.031 


64. 


.339 


600 


320 


0. 


15106 


16. 


.070 


30 


.047 


35 


.362 


63 


,953 


600 


330 


0. 


15042 


16 


.072 


30 


.118 


36 


.708 


63 


.752 


600 


340 


0. 


14968 


16 


.061 


30 


.168 


38. 


.071 


63 


,287 


600 


350 


0. 


15232 


16 


.340 


30 


.965 


39 


.528 


69 


,329 


600 


360 


0. 


14849 


16 


.053 


30 


.297 


40. 


.808 


62 


,629 


800 


180 


0. 


16884 


16 


.203 


29. 


.105 


12. 


.929 


70. 


.077 


800 


190 


0. 


16688 


16. 


.200 


29 


137 


13 


.823 


69 


.926 


800 


200 


0. 


16491 


16. 


.183 


29 


.142 


14. 


.726 


69. 


.348 


800 


210 


0. 


16455 


16 


.283 


29 


.351 


15. 


.607 


71. 


,847 


800 


220 


0. 


16132 


16. 


135 


29 


.140 


16 


.558 


67. 


,801 


800 


230 


0. 


15982 


16. 


.127 


29. 


.147 


17. 


.478 


67 


.534 


800 


240 


0. 


15858 


16. 


.118 


29 


.174 


18. 


.397 


67. 


.030 


800 


250 


0. 


15734 


16. 


.107 


29 


.193 


19 


.327 


66 


.550 


800 


260 


0. 


15624 


16 


.099 


29 


.217 


20. 


.258 


66 


.079 


800 


270 


0. 


15527 


16 


.092 


29. 


.249 


21, 


.193 


65 


.645 


800 


280 


0. 


15432 


16. 


.085 


29 


.276 


22 


133 


65 


,237 


800 


290 


0. 


15341 


16, 


.076 


29. 


.305 


23 


.083 


64 


.764 


800 


300 


0. 


15265 


16 


.073 


29. 


.342 


24. 


.028 


64. 


.400 


800 


310 


0. 


15186 


16. 


.067 


29 


373 


24. 


.983 


64. 


.001 


800 


320 


0. 


15117 


16 


.063 


29 


.414 


25 


.943 


63 


.644 


800 


330 


0. 


15047 


16 


.057 


29 


.448 


26 


.909 


63 


.268 


800 


340 


0. 


14983 


16 


.053 


29 


.489 


27 


.881 


62 


,933 


800 


350 


0. 


14925 


16 


.049 


29. 


.529 


28. 


.852 


62 


.575 


800 


360 


0. 


14866 


16. 


.041 


29 


.572 


29. 


.839 


62 


.127 
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Table A.l — continued 



-y 


K 




n 


- 


w 


Sq 


L 


^0 


1000 


180 


0. 


16891 


16. 


.213 


28. 


.929 


10.276 


70 


.446 


1000 


190 


0. 


16656 


16. 


.179 


28. 


890 


10.985 


69 


,460 


1000 


200 


0. 


16460 


16. 


.164 


28. 


.888 


11.700 


68. 


.969 


1000 


210 


0. 


16411 


16. 


.266 


29, 


106 


12.415 


71 


,720 


1000 


220 


0. 


16127 


16. 


.133 


28. 


.888 


13.136 


67 


.828 


1000 


230 


0. 


15989 


16. 


.124 


28. 


890 


13.853 


67 


.378 


1000 


240 


0. 


15857 


16. 


.114 


28. 


.907 


14.584 


66 


.911 


1000 


250 


0. 


15738 


16. 


103 


28. 


920 


15.314 


66 


.392 


1000 


260 


0. 


15626 


16 


.095 


28. 


.933 


16.047 


65 


,958 


1000 


270 


0. 


15525 


16. 


.086 


28. 


.950 


16.783 


65 


.496 


1000 


280 


0. 


15433 


16. 


.080 


28. 


969 


17.520 


65 


.071 


1000 


290 


0. 


15345 


16 


.074 


28. 


.988 


18.261 


64 


,677 


1000 


300 


0. 


15268 


16. 


.069 


29. 


.012 


19.001 


64 


.238 


1000 


310 


0. 


15191 


16 


063 


29. 


.037 


19.752 


63 


,864 


1000 


320 


0. 


15121 


16. 


.059 


29. 


060 


20.499 


63 


.503 


1000 


330 


0. 


15053 


16 


.052 


29. 


.085 


21.253 


63 


.069 


1000 


340 


0. 


14991 


16. 


.048 


29. 


110 


22.008 


62 


.729 


1000 


350 


0. 


14935 


16. 


.045 


29. 


.141 


22.763 


62 


.395 


1000 


360 


0. 


14880 


16 


.042 


29. 


.171 


23.525 


62 


,056 


1400 


180 


0. 


16872 


16. 


.195 


28. 


.677 


7.2841 


69 


.908 


1400 


190 


0. 


16657 


16 


.178 


28. 


.657 


7.7828 


69 


,419 


1400 


200 


0. 


16460 


16. 


.160 


28. 


.641 


8.2861 


68 


,842 


1400 


210 


0. 


16291 


16. 


.146 


28. 


633 


8.7877 


68. 


.295 


1400 


220 


0. 


16130 


16. 


132 


28. 


624 


9.2955 


67 


,785 


1400 


230 


0. 


15979 


16. 


.115 


28. 


.608 


9.8043 


67 


,148 


1400 


240 


0. 


15867 


16. 


.117 


28. 


633 


10.311 


66 


,953 


1400 


250 


0. 


15735 


16. 


.097 


28. 


.615 


10.825 


66 


.220 


1400 


260 


0. 


15632 


16. 


093 


28. 


.625 


11.336 


65 


.868 


1400 


270 


0. 


15527 


16. 


.084 


28. 


.627 


11.852 


65 


.403 
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Table A.l — continued 



—y 


Kn 

U 








Wq 


So 


L 


^0 


1400 


280 


0. 


15433 


16 


.075 


28. 


,634 


12.369 


64. 


.909 


1400 


290 


0. 


15352 


16 


.071 


28 


.645 


12.883 


64. 


.546 


1400 


300 


0. 


15269 


16 


,063 


28 


.655 


13.405 


64. 


.064 


1400 


310 


0. 


15196 


16 


.058 


28 


.667 


13.924 


63. 


663 


1400 


320 


0. 


15126 


16 


.053 


28 


.681 


14.447 


63. 


317 


1400 


330 


0. 


15065 


16 


.048 


28 


.693 


14.965 


62. 


.884 


1400 


340 


0. 


15003 


16 


,044 


28. 


.708 


15.490 


62. 


.512 


1400 


350 


0. 


14941 


16 


.039 


28 


.713 


16.015 


62. 


.182 


1400 


360 


0. 


14889 


16 


.035 


28 


.738 


16.544 


61. 


.778 


1800 


180 


0. 


16868 


16 


,193 


28 


.552 


5.6423 


69. 


.913 


1800 


190 


0. 


16652 


16 


,174 


28 


.525 


6.0275 


69. 


339 


1800 


200 


0. 


16460 


16 


,158 


28 


.508 


6.4145 


68. 


.779 


1800 


210 


0. 


16281 


16 


.142 


28 


,488 


6.8046 


68. 


.254 


1800 


220 


0. 


16126 


16 


.129 


28 


.475 


7.1940 


67. 


.714 


1800 


230 


0. 


15983 


16 


.116 


28 


,465 


7.5857 


67. 


193 


1800 


240 


0. 


15846 


16 


.104 


28 


.446 


7.9788 


66. 


.736 


1800 


250 


0. 


15738 


16 


.098 


28 


.459 


8.3715 


66. 


.256 


1800 


260 


0. 


15629 


16 


.090 


28 


,457 


8.7665 


65. 


790 


1800 


270 


0. 


15530 


16 


,082 


28 


,457 


9.1619 


65. 


.332 


1800 


280 


0. 


15437 


16 


.074 


28 


.459 


9.5591 


64. 


.869 


1800 


290 


0. 


15349 


16 


.065 


28. 


.457 


9.9567 


64. 


391 


1800 


300 


0. 


15271 


16 


,058 


28 


,462 


10.354 


63. 


.912 


1800 


310 


0. 


15199 


16 


.056 


28 


.471 


10.754 


63. 


.584 


1800 


320 


0. 


15129 


16 


,050 


28 


.478 


11.155 


63. 


.167 


1800 


330 


0. 


15068 


16 


.046 


28 


,485 


11.553 


62. 


.780 


1800 


340 


0. 


15009 


16. 


.043 


28 


.498 


11.955 


62. 


.431 


1800 


350 


0. 


14948 


16 


.037 


28 


.507 


12.361 


62. 


.062 


1800 


360 


0. 


14895 


16 


.032 


28 


.515 


12.763 


61. 


.648 
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